LAKE AND SHORE ICE CONDITIONS ON SOUTHEASTERN LAKE MICHIGAN IN THE 
VICINITY OF THE DONALD C. COOK NUCLEAR PLANT: WINTER 1973-74 



By 



Erwin Seibel, Christopher T. Carlson 
and Joseph W. Maresca, Jr.-'- 



Under contract with 

American Electric Power Service Corporation 

Indiana and Michigan Electric Company 

John C. Ayers, Project Director 



Special Report No. 55 
Great Lakes Research Division 
The University of Michigan 
Ann Arbor, Michigan 
1975 

'•Present address: Stanford Research Institute, Menlo Park, Calif. 94025 



TABLE OF CONTENTS 

ABSTRACT v 

ACKNOWLEDGMENTS vi 

INTRODUCTION 1 

METHODOLOGY 3 

GEOMETRY OF OBLIQUE PHOTOGRAPHS 7 

DISCUSSION AND RESULTS 10 

Terminology 10 

Ice Conditions and Climatological Data 11 

Ridge Location and Nearshore Breaker Zones 24 

Nearshore Erosion 28 

CONCLUSION 28 

REFERENCES 29 

APPENDIX A. DAILY ICE CONDITIONS 31 

APPENDIX B. DESCRIPTION OF SLIDES USED IN ANALYSIS 37 

APPENDIX C. SOURCES OF POSSIBLE ERRORS IN THE PHOTOGRAPHIC 

ANALYSIS 43 

APPENDIX D. METHOD USED IN THE DETERMINATION OF 

MONITORED VARIABLE POSITIONS 48 

APPENDIX E. COMPUTER PROGRAM ICESTUBI 50 

APPENDIX F. SUMMARY TABLE OF CLIMATIC VARIABLES: 

JANUARY AND FEBRUARY 1974 59 



ABSTRACT 

A time lapse photographic system designed to take high oblique photographs 
of the nearshore zone at the Donald C. Cook Nuclear Plant in southeastern Lake 
Michigan was used to provide a nearly continuous record of ice conditions 
during the 1973-74 winter season. The ice conditions have been categorized 
into five distinct stages: no ice, static, accretion, deterioration and break- 
up. These five stages of ice development were related to wind direction, wind 
speed, air temperature and water temperature. A typical sequence of shore ice 
formation through breakup was established through qualitative analysis of the 
photographs. The typical sequence was influenced by the complex interrelation- 
ship between the climatic variables and ice stage development. Quantitative 
analysis of the photographs was used to test and subsequently verify the hy- 
pothesis that the nearshore ice ridges, offshore bars, and breaker zones are 
coincident in location at approximately 40 and 100 meters from the water's edge. 
The observations reveal that large quantities of sediment are incorporated into 
the nearshore ice and that the nearshore ice ridges ground themselves on the 
nearshore bottom. 

We believe that the grounded nearshore ice ridges simultaneously modify 
the nearshore topography and protect the shoreline and bluffline from erosion 
by winter storms. The degree to which the nearshore is protected or modified 
has not yet been established. 
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INTRODUCTION 

As part of a continuing environmental study of southeastern Lake Michigan 
in the vicinity of the Donald C. Cook Nuclear Plant (Fig. 1), a series of time 
lapse photographs were taken during the winter of 1973-1974. The high oblique 
photographs depicting the shore and nearshore zones north of the plant provided 
a nearly continuous record of ice conditions during that period (see Appendix 
A) . This made possible the observation of shore ice formation and breakup as 
a function of climatological variables and nearshore dynamics. 

Due to the short term or irregular nature of their observations, early 
observers of shore ice characteristics were usually restricted to qualitative 
investigations of morphology or process. Some investigations (Heap 1963; Heap 
and Noble 1966) were dependent upon the expertise of various shore based ob- 
servers at numerous stations over a large geographical area, while others 
(Marshall 1966) based their studies on aerial photo reconnaissance. Such data 
gathered for the gross evaluation of the Great Lakes ice distributions and 
durations of coverage could not be used in the determination of smaller scale 
localized formative processes. 

More recently attention has been focused on the relationship between 
nearshore ice and nearshore morphology (Marsh et al. 1973; Ayers and Yocum 
1972) and the genesis of shore ice in response to various climatological, 
geomorphological and hydraulic parameters (O'Hara and Ayers 1972; Bryan and 
Marcus 1972) . Evenson (1973) reviewed efforts by Clos-Arceduc and Sonu to 
model the apparent regularity of beach cusps. Evenson suggested that if a 
model such as that developed by Clos-Arceduc was operational in the Great Lakes 
it might also explain the observed regularity of blowhole cones which are fre- 
quently aligned along the lakeward edge of the ice ridge complex. Attempting 
to establish a sequence of processes typical of shore ice formation in their 
geographical area of study, O'Hara and Ayers (1972) based their conclusions on 
the observations of photographs supplemented by intermittent field observa- 
tions. Ayers et al. (1973) initiated a time lapse photographic monitoring 
system into their study program to supplement aerial photography. Marsh et al. 
(1973) utilized seismic refraction and drilling to correlate nearshore ice 
structure with subsurface morphology. 



MILWAUKEE 




CHICAGO 



Ktlomatare 



JOSEPH 
/COOK PLANT 



'MICHIGAN CITY 



FIG. 1. Location of the Donald C. Cook Nuclear Plant. 



O'Hara and Ayers (1972), Ayers et al. (1973) and Bryan and Marcus (1972) 
proposed that there was a correlation between the location of ice ridges and 
offshore sand bars. O'Hara and Ayers (1972) obtained evidence in support of 
this supposition through diving operations conducted on 10 March 1970 during 
which grounded remnant ice blocks of a deteriorated ice ridge were observed 
to be embedded in the bottom and thought to sit on one of the offshore bars. 
Ayers et al. (1973) superimposed photographs of breaking waves and ice ridges 
which provided additional qualitative evidence of the ice ridge - offshore bar 
correlation. During an aerial reconnaissance flight Bryan and Marcus (1972) 
observed isolated ice blocks, remnants of a def -iorating ice ridge, grounded 
along a continuous offshore bar. 

Offshore profiles compiled by Marsh et al. (1973) revealed that a cor- 
relation does exist between offshore bars and ice ridges but they observed 
that within 150 meters from the Lake Superior shoreline bars existed without 
ice ridges located over them. 

In this study, the writers investigated the correlation between breaker 
zone position, offshore topography and ice ridge location. Computer analysis 
of time lapse high oblique photographs provided a quantitative basis for dis- 
cussion. Meteorological data were utilized in conjunction with the photographs 
to evaluate ice ridge development as a function of time, air and water temper- 
atures and wind speed and direction. 



METHODOLOGY 

High oblique photographs were obtained for analysis with a Canon F-1 
single reflex 35 mm camera with Canon 50 mm f 1.8 lens. The camera with Canon 
EE Servo Finder and special timing unit was mounted indoors facing in a NNW 
direction and automatically provided an average of five photographs per day 
from early morning to late afternoon. Since it was not possible to determine 
the exact time of exposure for each slide, a relative time was assigned to 
each photo: early am, late am, noon, early pm, and late pm. The positions 
of a sheet pile seawall and a range pole, evident in each photograph, were 
determined by standard field survey methods and were used as reference points. 



The camera location permitted observation of isolated ice blocks as well as 
cross -sectional ice ridge profiles visible through wave breached portions of 
the ridge system. In each slide the horizon was evident. Reflected glare 
from the lake surface, ice, or window glass did not interfere with the quality 
of the photographs. 

Slides were selected for use in the analysis instead of prints for two 
reasons. First, slides are less expensive to process than prints; and 
secondly, analysis required that the images of variables in the oblique 
photographs be marked on a tracing paper overlay. Slides were easily enlarged 
with a slide projector to minimize human error in selecting points from the 
image. 

Representative slides or groups of slides were chosen for analysis from 
the sequence taken during the period from 17 January 1974 through 6 April 1974 
(Appendix B) . Although the camera was set up late in December 1973, photo- 
graphs were not available until 17 January 1974 due to camera malfunction. As 
a consequence, the initial ice formation occurring late in December 1973 was 
not recorded. The final melt-off of Isolated remnant ice blocks occurred on 
6 April 1974. 

Analysis of the slides was divided into three steps. First, points were 
marked locating the positions of variable features in the projected oblique 
representation, and their coordinates were determined relative to an (x, y) 
coordinate system with origin at the principal point of the oblique photo- 
graph. Secondly, the positions of these points were calculated for the 
equivalent vertical photograph with coordinate system origin at its isocenter. 
Finally, the real ground coordinate positions of the points were determined 
relative to a coordinate system with origin at the camera. The geometry basic 
to this method of oblique photo analysis is discussed under Geometvy of 
ohtique photographs. 

A fixed reference system was requisite for slide-to-slide comparisons. 
Establishment of such a coordinate system with origin at the principal point 
of the projected oblique photograph necessitated the removal of the color 
transparency from the developer's cardboard mount and remounting on a 2,0 inch 
plexiglass square as suggested by Maresca (1975). Inscribed on the plastic 
mount were a series of lines parallel to its borders to aid in placement of 



the positive. Also etched on the mount were the fiducial axes whose origin 
would be located at the principal point of the remounted color transparency. 

The remounted slide was projected by a Kodak Ektagraphic slide projector 
onto a wall covered with a tracing paper overlay. The projected image of the 
fiducial coordinate axes, now superimposed on the projected photograph, was 
located on the overlay. To minimize error due to remounting, the location of 
the projected principal point was checked and adjusted when necessary. 

The perpendicular distance was measured between two of the inscribed 
horizontal lines on the plastic mount, one at the top and one at the bottom 
of the mount. The distance between the projected image of these same two 
lines was measured on the screen, and these values were used in the computer 
program for calculation of the focal length of the projector. The positions 
of the reference points common to all slides were marked on the overlay and 
subsequently the projected positions of the variables unique in each slide 
were located. These variables included the horizon, breaker zones, ice ridges 
and interridge ice lagoons, zones of brash ice accumulations, heights of the 
icefoot, ice blocks, and ice ridges and the waterline when visible. 

The computer program utilizes the mean still water level (MSWL) as the 
datum for ground coordinate calculations thus requiring that the selected data 
points be located at this elevation. Points defining variable positions were 
therefore marked on the tracing paper overlay where their intersection with 
the water surface was evident in the projected picture. The points marking the 
location of breaker zones were placed at the base of the breaking waves. 
Points marked along the crest of an ice ridge, above the MSWL datum, would 
represent only an estimate of the ridge's position (see Appendix C) . Conse- 
quently, when determining ice ridge location, it was frequently necessary to 
use a later slide which depicted the breaching of or breakup of that ice ridge 
when its intersection with the water surface at midridge was visible. Points 
defining the outline of brash ice areas of the edges of the interridge ice 
lagoons were in most cases considered to be at the MS'WL. The calculation of 
the height of a variable required two points of definition: one at the 
object's base where intersection with the water was apparent and the other 
marking its vertical extent. For the reference points, which were not located 
at the still water level, surveyed heights were specified in the program for 



determination of their ground coordinate location. 

The distance from the principal point of the oblique photograph to the 
apparent horizon is utilized in the program for calculation of distances on the 
principal plane diagram, hence, the position of a poorly defined horizon may 
introduce errors. The quality of the horizon was evaluated as excellent, good, 
fair, or poor, thus offering one possible explanation of error in poor computer 
output values. 

The determination of the (x, y) coordinates for the points locating the 
positions of variables on the projected oblique photograph required a coordinate 
grid. Coordinate values were read to the third decimal place from a cartesian 
coordinate grid with a scale of twenty divisions to the inch placed under the 
tracing paper overlay. 

The computer program converts the oblique (x, y) coordinate points to co- 
ordinates for an equivalent vertical photograph and then calculates the real 
ground coordinates for each ordered pair. Printer plot subroutines, called in 
the main program, construct a ground coordinate grid with origin at the camera. 
The computed ground coordinates of each data point were scaled and plotted 
within the grid. A unique plotted character was used for each monitored vari- 
able. Characters representing the location of each variable were positioned in 
linear arrays on the printer plots and were approximated by straight lines to 
facilitate distance measurements as described in Appendix D. The distances to 
all variable positions were calculated perpendicular to a baseline along the 
shore. The sheet pile seawall was chosen as the baseline since the water line 
is subject to frequent changes in position and is usually not visible during 
periods of ice cover. Presently the computer program does not calculate the 
distances from the baseline to the various variable positions. A reference 
point along the baseline was established and all distance measurements were 
made perpendicular to the baseline from that point. A normal to the baseline 
was constructed of sufficient length to intersect the lines defining the vari- 
able positions. This normal formed the hypothenuse of a right triangle. 
Determination of the x and y components of the baseline normal and application 
of the plotting scale factor and the Pythagorean Theorem yielded distances to 
the variable's location. This technique was used to calculate the distances 
from the baseline to all monitored variables. 



GEOMETRY OF OBLIQUE PHOTOGRAPHS 

The analysis of an oblique representation can be obtained in any standard 
reference or text (American Society of Photo grairanetry 1966) , but is presented 
here briefly so that the correction for the projection of the image through a 
slide projector can be included. The geometry of the principal plane diagram 
(Olson 1973) (Fig. 2) was used to derive the fundamental equations required to 
solve the oblique photograph problem. The distance, PT, is measured on the 
screen from the apparent horizon to the principal point of the photograph. 
Since this distance is projected onto a screen from a slide projector, the 
image distance must be corrected for the focal length of the slide projector. 
A constant, C, is used to make this correction. The constant, C, is defined as 

C = I/O (1) 

where I is the image distance and is the object distance. The image 
distance, I, is calculated from the lens maker's formula 

I = Oxf./0-f. (2) 

pro3 proj 

where f . is the focal length of the projector. The distances on the 
principal plane positive can now be defined. The angle of depression, 9', is 
defined as 



tatx 9' = PT / f (3) 

cam 

where PT = PT x C (4) 

screen 

and f is the focal length of the camera. The distance PT and the angle 9' 
cam * 

was calculated to the apparent horizon. A small correction can be made to 
determine 9' to the true horizon by adding 

AS = 0.98 /¥ (5) 
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FIG. 2. Principal plane diagram. 



in minutes to 6' where H is the difference in elevation between the water 
surface and the camera elevation. Since 9 = 6 ' + A9 , the distance PT must be 
re-calculated. The geometry of the principal plane diagram (Fig. 2) may be 
used to derive the distance PI using 

PI = (tan t/2) X f (6) 

cam 

where t equals 90-9. The remaining distances in the principal plane diagram 
are given in equations 7, 8, and 9. 

PN = f X (tan t) (7) 

cam 

TI = PT + PI (8) 

TN = PT + PN (9) 

Using this information all points in the oblique representation can be rotated 

into an equivalent vertical photograph and finally real ground distances were 

obtained. 

Individual ordered pairs, (X , Y ), are picked off the oblique photo- 

scrn scrn ' '^ -i r 

graphic projection on the screen relative to a rectangular coordinate system 
with an origin at the principal point. The coordinates, (X , Y ) of the 
equivalent vertical photographs are calculated as 

X = X X C X (TI / (TI - K)) (10) 

v.p. scrn 

Y = K X (TI / (TI - K)) (11) 

P * 

where K = (Y x C) + PI (12) 

scrn ' 

The real ground coordinates, (X j» Y ), can be calculated from a simple pro- 
portion between the equivalent vertical photograph and the ground relative to 

a rectangular coordinate with the origin at the camera. The (X ,, Y ,) is 

^ grd' grd 

calculated as 



X,=HxX /f (13) 

grd v.p. cam ^ ' 

Y ,=Hx(PI+Y )/f (14) 

grd ^ v.p. cam ^ '' 

The fundamental assumption basic to the analysis is that all points lie within 
the datum plane. If a point lies outside the datum, as for example a point on 
the beach, then a different elevation, H, must be used. 

Vertical elevations or vertical heights were calculated from the formula 

d' TN 
h = H X -^ X ^ (15) 

where d' is the height of the object on the screen, r' is the distance from the 
principal point to the top of the object in the oblique view, TN is the distance 
from point T to point N, and TS is the distance from the point T to the point S 
which is the base of the object. 

Using the above formulae, a comprehensive computer program, ICESTUDY 
(Appendix E) , was written to convert from the oblique to the real ground 
distances. The output was produced in both tabular and graphical format. The 
breaker zone and ice ridge locations were calculated as well as the position of 
the wave runup. The height of the ice ridges was calculated when visible in 
the photo and the area of zones of accumulated brash ice was determined . 

DISCUSSION AND RESULTS 

TERMINOLOGY 

The terminology used in this report is that of O'Hara and Ayers (1972) and 
is consistent with that recommended by Klvisild (1970). The initial stage of 
lake ice development is the formation of a frozen beach. This is usually 
followed by the buildup of a small ridge or icefoot at or near the water's edge 
(Fig. 3a). The icefoot is the result of spray produced by waves lapping against 
the frozen beach face and the accumulation of brash ice fragments. Subsequent 
to the icefoot development is the formation of a tabular ice lagoon which ter- 
minates at the first ice ridge (Fig. 3b). Numerous, small, irregularly spaced 
undulations may be observed on the ice lagoon and are oriented roughly parallel 
with the lakeward edge of the accreting ice mass. Each undulation marks a 
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short hiatus in the lakeward ice advance when freezing spray and minor amounts 
of wave tossed brash ice accumulate on the advancing icefront. The first ice 
ridge is generally the smallest of the lake ice ridges. Like its more lake- 
ward counterparts (Fig. 3c) it has a steep nearly vertical lakeward face and a 
gently sloping landward surface. Both ridges are frequently reduced to 
grounded remnant ice blocks upon deterioration or bre£ikup (Fig. 3d). For the 
discussion that follows, the ice complex at the study site is defined by the 
patterned combination of three features: the icefoot, the ice lagoon and the 
ice ridge. A schematic representation of these three structures is shown in 
Figure 4. 

ICE CONDITIONS AND CLIMATOLOGICAL DATA 

The ice conditions were categorized into five distinct stages: no ice, 
accretion, static, deterioration and break up. The stages of ice development 
used here were determined from the time lapse photography of the ice at the 
study site. Accretion is a visible addition of ice since the time of the 
previous slide while deterioration and breakup are a reduction in the ice mass. 
Deterioration is a gradual reduction of the ice mass apparent only after a 
series of slides is viewed in succession while breakup is the rapid destruction 
of the nearshore ice complex usually under severe lake conditions. During 
breakup, rapid change is visible from slide to slide and is usually exemplified 
by the landward displacement of the ice complex. The static condition referred 
to here connotates no visible change in the ice complex from slide to slide 
while the condition of "none" refers to the total absence of ice in the picture. 

The stages of ice development were related to the following climatic con- 
ditions: wind direction, wind speed, air temperature and water temperature. 
Figures 5 through 10 and Appendix F show the climatic conditions as they relate 
to the ice conditions for January and February 1974, There does not appear to 
be a single climatic variable that can be correlated alone to ice development 
stages. The combined influence of air temperature, wind direction and wind 
speed control the stages of ice development in a complicated manner. 

It may be considered an obvious observation that increasing or decreasing 
the air temperature directly influences the breakup or accretion of lake ice. 
During the months of January and February 1974 ice deteriorated or began to 
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FIG. 3. Characteristic ice features at various stages in the for- 
mation and breakup of the near shore ice complex. 
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c) The ice complex from shore to beyond the 
second ice ridge. 
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d) Remnant ice blocks delineating the former 
position of the second ice ridge. 

FIG. 3 continued. 
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breakup shortly after an increase in air temperature. Subfreezing air temper- 
atures usually initiated accretion; however, during the period 14-15 February 
1974 ice breakup occurred despite subzero air temperatures. During this time 
period the water temperature remained constant but the wind direction was pre- 
dominantly offshore at speeds of about 10 knots, indicating that the accre- 
tionary effect of subzero air temperatures may be negated by strong offshore 
winds . 

The climatic variable that has the least influence on the stages of ice 
development defined here is the surface water temperature. Figure 6b and 
Figure 9b represent the surface water temperatures and the ice development 
stages for January and February 1974 respectively. An examination of these 
figures indicates a poor correlation between these variables. 

To more closely examine the relationship between the climatic variables 
and the ice conditions we chose an accretionary sequence starting on 22 
February and extending through 26 February, for which the time lapse system 
provided good coverage. For this period Figure lla-d shows the ice conditions 
plotted against each of the climatic variables, and Figure 12 represents the 
ice buildup from the baseline. Several observations can be made that may lead 
to a better understanding of the nearshore ice development at the study site. 

Immediately preceding the accretionary sequence, easterly winds and rising 
air temperatures produced significant overnight deterioration on 21 February 
1974 and reduced the preexisting ice ridge system to grounded blocks and 
floating ice fragments. Rotation of wind direction to an onshore northwesterly 
quadrant accompanied by an increase in wind speed and a drop in air tempera- 
tures on 22 February 1974 initiated the accretive sequence. The onshore wind 
direction coupled with high wind speed forced the remnant ice blocks onto the 
beach and subsequently augmented this accumulation by moving floating brash, 
ball and pancake ice shoreward, producing a nearly vertical 2 meter high wall 
of shore ice in less than 3 hours. Rapid ice accumulation under severe lake 
conditions on 22 February resulted in the incorporation of large volumes of 
sediment in the ice giving it a gray-brown coloration. Ice formed under the 
less severe conditions during the remainder of the accretive sequence (23-26 
February) was much whiter in color, making it easily distinguishable from the 
darker colored ice deposited during the storm. Although air temperatures were 
well below freezing. Figure 11 shows one data point (24 February) indicative 
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of significant deterioration. This corresponds to a photograph taken just 
following a period of approximately 8 hours when wind direction was pre- 
dominantly offshore, again demonstrating the complex interrelationship between 
ice stage development and climatic variables. 

RIDGE LOCATION AND NEARSHORE BREAKER ZONES 

One of the primary aspects of our work was oriented to the quantification 
of the ice ridge position during the winter season and its location relative to 
the nearshore breaker zones. Ayers et al. 1973 suggest that visual observation 
and the qualitative superimposition of slides depicting the location of the 
breaker zones with that of the ice ridges leads to the conclusion that the two 
are coincidental. 

Using the oblique photographs for the winter of 1974, we measured the dis- 
tance from a permanent baseline to the breaker zones and the ice ridges. The 
measured points are shown in Figure 13 and their values recorded in Table 1. 
The coincidence of the location of the breaker zones and the ice ridges is 
apparent. Three breaker zones and three ice ridges were identified during 
this winter period. The ice ridge nearest to shore and its corresponding 
breaker zone are best defined with the least scatter in the data points. As 
the distance from shore increases, coincidence between the breaker zone and 
ice ridges is still present but the scatter becomes much greater as summarized 
in Table 2. 

The scatter in the data is a function of the analysis and the real varia- 
bility in the mean position of the ice ridge, breaker zone and nearshore bars. 
Each variable is positioned within a linear zone and its location is approxi- 
mated by a line for measurement purposes. In particular, the breaker zone 
(offshore bar) has a natural width and since its location is defined by only 
one slide at a given time, this location is dependent upon wave conditions as 
well as the position of the breaking wave at the time of photo exposure. It 
follows that one photo at any given time is not sufficient to define the range 
in width of this zone. The position of an ice ridge is more definitive since 
grounded remnant ice blocks (Fig. 3d) maintain their position during breakup. 
Frequently wave activity breaches an ice ridge providing a cross-sectional 
view and revealing its mid-ridge intersection with the water surface. 
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photographs. 
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TABLE 1. Identification of the case numbers used in Figure 13. 









0-75 Meters 


76-150 ] 


Meters 


151-300 Meters 








Breaker 


Ice 


Breaker 


Ice 


Breaker Ice 


Case 


Slide No. 


Date 


Zone 


Ridge 


Zone 


Ridge 


Zone Ridge 


1 


CM-9-74-23 


1-17-24 








114 




2 


CM-10-74-2 


1-22-74 






114 






3 


CM-10-74-3 


1-22-74 




45 








4 


CM-10-74-9 


1-23-74 


41 




88 




172 


5 


CM-10-74-18 


1-24-74 


47 










6 


CM-10-74-20 


1-25-74 




49 


98 






7 


CM-10-74-22 


1-25-74 




43 


98 






8 


CM-10-74-2 3 


1-25-74 




47 








9 


CM-10-74-28 


1-26-74 




46 








10 


CM-10-74-29 


1-27-74 


46 




105 






11 


CM-11-74-5 


1-29-74 


48 










12 


CM-11-74-6 


1-30-74 


48 










13 


CM-11-74-13 


1-31-74 


45 




127 




215 


14 


CM-11-74-30 


2-3-74 


55 




105 






15 


CM-11-74-32 


2-4-74 


45 




89 






16 


CM-12-74-12 


2-7-74 


48 




94 






17 


CM-12-74-18 


2-8-74 




46 








18 


CM-12-74-20 


2-8-74 




46 








19 


CM-12-74-29 


2-10-74 




50 








20 


CM-13-74-13 


2-14-74 




45 








21 


CM-13-74-19 


2-15-74 










216 


22 


CM-13-74-24 


2-16-74 




53 








23 


CM-13-74-26 


2-17-74 








112 




24 


CM-14-74-2 


2-19-74 








110 


218 


25 


CM-14-74-3 


2-19-74 






113 


123 


237 


26 


CM-14-74-15 


2-22-74 




44 








27 


CM-14-74-16 


2-22-74 


50 




105 




177 


28 


CM-14-74-17 


2-22-74 


37 




109 






29 


CM-14-74-18 


2-22-74 


49 




100 




214 


30 


CM-14-74-19 


2-22-74 


47 




110 




176 


31 


CM-14-74-21 


2-23-74 


55 








153, 270 


32 


CM-14-74-22 


2-23-74 


62 




118 




194 


33 


CM-14-74-23 


2-23-74 


54 




128 






34 


CM-14-74-24 


2-23-74 


38, 59 




134 






35 


CM-14-74-26 


2-24-74 


55 










36 


CM-14-74-27 


2-24-74 


52 










37 


CM-14-74-28 


2-24-74 


57 




144 






38 


CM-14-74-2 9 


2-24-74 


59 




136 






39 


CM-14-74-30 


2-25-74 






116 






40 


CM-14-74-31 


2-25-74 


57 




100 






41 


CM-14-74-33 


2-25-74 






123 






42 


CM-14-74-34 


2-25-74 






109 






43 


CM-14-74-36 


2-26-74 








115 




44 


CM-17-74-3 


3-13-74 


51 




113 






45 


CM-18-74-5 


3-20-74 


46 
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TABLE 2. Statistical summary of breaker zone and ice ridge locations. 



y 


DISTANCE FROM BASELINE (METERS) 


y^ 


0-75 


76-150 


151-300 


y/^ 


n 


X 


Sx 


n 


X 


Sx 


n 


X 


Sx 


BREAKER 
ZONES 


25 


50 


7 


24 


112 


15 


10 


203 


35 


ICE 
RIDGES 


11 


47 


3 


5 


115 


5 


1 


216 





n = number of cases 



X = 
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i=l 


^i 
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o2 


n 
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i=l 


(x.2) - 


n 


2 
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b — 

X 






n - 


- 1 







This breached cross-section maintains its position in several slides allowing 
calculation of both the ice ridge's position and height. 

The position of the breaker zones is known to reflect the position of the 
offshore bar system. The clear relationship between the offshore ice ridges 
and the breaker zones leads to the conclusion that ice ridges are coincidental 
with the bar systems. The time lapse photographs of the nearshore ice complex 
indicate that once established the ice ridges remain stationary for longer 
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periods than the intervening lagoons. This suggests that not only are the ice 
ridges greater in thickness than the ice lagoons but that they are situated 
directly on the nearshore lake bottom. On one occasion, O'Hara and Ayers 
(1972) verified during diving operations that the remnant blocks of a deteri- 
orating ice ridge were embedded in the bottom. Bryan and Marcus (1972) ob- 
served remnant ice blocks grounded along offshore bars during an aerial recon- 
naissance. 

NEARSHORE EROSION 

Seibel (1972) formulated that the ice ridge formation and ice ridge 
breakup are an integral part of the modification of offshore topography and 
aggravation of shoreline erosion. Scour of the nearshore topography is in- 
evitable at the position of the ice ridges. The influence of the modification 
on the shoreline and bluffline is unknown. It is speculated that a storm 
immediately following breakup could cause greater damage along the shoreline 
than a storm of equal magnitude several weeks after breakup when the nearshore 
bars have had a sufficient opportunity to reform. 

It has been clearly observed by ourselves and others that large quantities 
of sediment are incorporated in the nearshore ice complex during accretionary 
stages. At the time of breakup this ice-locked sediment was transported land- 
ward, lakeward or alongshore. The combined effect of ice rafting and scour 
suggest that nearshore ice plays a greater than anticipated role in modifying 
the nearshore topography. 



CONCLUSION 

Several pertinent conclusions appear to be evident from the monitoring of 
the nearshore ice on this portion of southeastern Lake Michigan: 

1. The computer analysis of the time lapse high oblique photographs pro- 
vided the basis for a quantitative analysis of the hypothesis of coincident 
location of offshore ice ridges, offshore bars, and breaker zones. 

2. The data reveal that the three above characteristic shore features 
are indeed coincident. 

3. The stages of ice development are not controlled by any single 
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climatic variable but instead there exists a complex interrelationship between 
the two. 

4. Easterly winds coupled with rising air temperatures appear to con- 
sistently produce significant deterioration of the ice complex. 

5. Onshore northwesterly winds coupled with increase in wind speed and 
a drop in air temperature initiate accretive sequences. 

6. Large quantities of sediment are incorporated in the nearshore ice 
complex during accretionary stages and are transported landward, lakeward or 
alongshore at the time of breakup. 

7. Since the ice ridges are grounded, as the data indicate, they may 
serve to modify the nearshore topography in the proximity of the offshore bars 
sufficiently to influence the wave regime capable of reaching and acting on 
the shoreline and bluff line. The exact influence of this modification of the 
shoreline and bluff line is not yet established. 
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APPENDIX A. DAILY ICE CONDITIONS 



17 January 1974 
(CM-9-74-23) 



The first photo taken from the came'.ra monitor system and 
utilized in this study was exposed 17 January 1974, Evident 
in the photo was the deteriorated first ice ridge, dis- 
colored due to sedimentary lag deposits. The second lagoon 
had been reduced to floating brash ice and the position of 
the second ridge was marked by grounded remnant ice blocks 
of average height 1-4 m. The remnant ice fragments of the 
second lagoon are absent in photos taken later on 17 January 
1974. 



18 January 1974 



There was continued deterioration of the first ridge and 
the remnant blocks of the second ice ridge. 



19-21 January 1974 No photos were available for these three days. 



22 January 1974 
(CM-10-74-2,3) 



23 January 1974 
(CM-10-74-9) 



24 January 1974 
(CM-lD-74-18) 



25 January 1974 
(CM-10-74-20,22, 
23) 



26 January 1974 
(CM-10-74-28) 



27 January 1974 
(CM-10-74-29) 



Waves displaced the ice blocks fornierly marking the position 
of the second ridge and rapidly reduced the first ridge to 
a narrow band of intermittently breiached, groimded ice. 

Overnight wave activity further reduced the first ridge to 
ice fragments and ice blocks of average height 1.3 m 
grounded along the shoreline. Wave action and continued 
influx of brash ice rebuilt the ice foot at the water's 
edge. 

Little change in the ice complex since 23 January. The 
lakeward edge of the icefoot had become more irregular in 
plan and the few isolated blocks of ice formerly located 
between the plunge zone and the icefoot were no longer 
present. 

Overnight ice accretion resulted in the formation of the 
first lagoon and the first ice ridge. The lagoon is com- 
posed of unconsolidated brash ice isolated behind the wave 
resistant first ridge. A band of slush ice was present 
lakeward of the first ice ridge. The first ridge was 
augmented by freezing spray and wave tossed slush ice which 
caused it to narrow and increase in height. 

The band of slush ice lakeward of the first ice ridge dis- 
sipated overnight and the deterioration of the lakeward 
edge of the first ridge accented its irregularity in plan. 
The icefoot visibly deteriorated. 

Overnight wave action broke up the first ridge and piled the 
ice fragments on the shore rejuvenating the icefoot. The 
intermittent patches of slush ice present in photos taken 
early 27 January 1974 were absent in later photos taken 
that day and the icefoot showed visible wave induced deteri- 
oration as the day progressed. 
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28 January 1974 



The ice conditions remained visibly unchanged other than a 
light snow cover. 



29 January 1974 
(CM-11-74-5) 



Wave action eroded the lakeward edge of the icefoot back to 
a position behind the waterline. This deterioration con- 
tinued until both the waterline and a narrow, irregular 
strip of beach were visible behind the icefoot in the last 
photo taken 29 January 1974. 



30 January 1974 
(CM-11-74-6) 



Approximately 21 m of overnight fast ice accretion advanced 
the lakeward edge of the fast ice to a location just shore- 
ward of the first breaker zone, no ice ridge was present. 
There was little further change in the ice conditions on 30 
January 1974. 



31 January 1974 
(CM-11-74-13) 



1 February 1974 



2 February 1974 



Northwesterly winds broke-up the fast ice overnight re- 
building the icefoot, the only ice feature visible in the 
first photo taken 31 January. The intense wave activity 
during the balance of the day finally reduced the icefoot 
to fragments which were washed onto the beach. 

The only remaining ice consisted of the remnant fragments 
of the icefoot covered by a light overnight snow. 

Light wave conditions and the constant influx of small 
quantities slush ice resulted in the slow formation of a 
small icefoot approximately 0.5 m high as the day progressed. 



3 February 1974 
(CM-11-74-30) 



The first two photos taken 3 February documented the con- 
tinued growth of the icefoot and the formation of small ice 
lobes on its lakeward edge. The increase in wave activity 
later in the day deteriorated the icefoot and produced a 
very irregular shape in plan. The icefoot was discolored 
due to the incorporation of sediment in the ice. 



4 February 1974 
(CM-11-74-32) 



Approximately 18 m of overnight fast ice accretion but no 
ridges were identifiable in early photos taken on this date. 
The height of the lakeward edge of the fast ice was in- 
creased by wave tossed brash ice, freezing spray generated 
by waves striking the ice formation and from recently de- 
veloped blowholes. 



5 February 1974 

6 February 1974 



7 February 1974 
(CM-12-74-12) 



No visible change in the ice conditions during this day. 

Although wave Intensity was low, wave splash and wave tossed 
slush ice slowly increased the height of the lakeward edge 
of the fast ice. New ice along the edge was darker than 
the rest of the snow covered fast ice indicating sediment 
incorporation into the ice. 

Approximately 8m of fast ice accretion overnight . The 
newly accreted ice was much whiter in color than the former 
lakeward edge of the fast ice making the boundary easily 
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8 February 1974 
(CM-12-74-18,20) 



9 February 1974 



10 February 1974 
(CM-12-74-29) 



11 February 1974 



12 February 1974 



distinguishable. Numerous blowhole cones were visible 
in photos taken on this date. 

Slush ice was visible near the horizon in the first photo 
taken 8 February and brash ice had begun to accumulate in 
a narrow band approximately 19 m wide just lakeward of the 
fast ice, but separated from it b3^ a zone of open water. 
The last photo taken 8 February showed that slush ice 
covered the entire lake and the band of accumulating brash 
ice had widened and consolidated. The narrow band of semi- 
consolidated brash ice remained isolated from the fast ice 
by the intermediate zone of the slush ice - water mixture. 

The entire lake was ice covered in the first photo taken 
this day. By midday drifting pancake ice was produced by 
the breakup evident in a zone lakeward of the main body of 
fast ice. The last photo taken 9 February showed the zone 
of floating brash ice had begun to reconsolidate. 

The lake was ice covered to the horizon in photos exposed 
early in the day, later photos show the horizon was icefree. 
A small ice ridge had begun to form over the former position 
of the first breaker zone, the location previously occupied 
by the band of brash ice evident in the photos taken 8 
February. 

Open water was visible near the horizon in the first photo 
taken on 11 February although later photos showed the lake 
was ice covered to the horizon. There was overnight enlarge- 
ment of the new ridge which had become evident in slides 
taken 10 February. 

A zone of open water formed near the horizon, its size did 
not appear to vary appreciably during the day. Otherwise 
there was no visible change in the other ice structures. 
Sun rot had continued to concentrate lag deposits around 
blow holes in the fast ice. 



13 February 1974 



14 February 1974 
(CM-13-74-13) 



The area of ice-free water evident near the horizon in 
photos taken 12 February was barely visible due to both its 
reduction in size and the foggy conditions which existed 
during the time of the midday photo exposure. A narrow 
strip of open water separated the fast ice from the shore 
at the water line. This may have been an atypical feature 
resulting from minor quantities of effluent discharge 
below the field of view at the bottom of the photos. 

The ice broke up into pancake and floating brash ice, re- 
sulting in a wide zone of open water bounded by floe ice 
near the horizon and the fast ice at the shore. Ice blocks 
were grounded along a line marking the position occupied by 
a breaker zone in previous slides. 
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15 February 1974 
(CM-13-74-19) 



16 February 1974 
(CM-13-74-24) 



17 February 1974 
(CM-13-74-26) 



18 February 1974 



19 February 1974 
(CM-14-74-2,3) 



20 February 1974 



21 February 1974 



22 February 1974 
(CM-14-74-15,16, 
17,18,19) 



The fast ice and grounded ice blocks were still present in 
photos taken this day. The lake was otherwise ice-free to 
the horizon since the floe ice located near the horizon in 
photos taken 14 February had been dispersed overnight by 
easterly winds. 

The first photo taken this day showed brash ice had begun 
to accumulate in a band adjacent to the lakeward edge of the 
fast ice. Waves forced the brash ice shoreward forming a 
small ridge located over the position of the first breaker 
zone. The newly forming ridge initially remained isolated 
from the fast ice by a zone of relatively ice free water but 
was later forced nearer to the fast ice by increased wave 
action. The grounded blocks of ice became more scattered by 
the end of the day and appeared to have moved in a northerly 
longshore direction. 

Overnight ice accretion of approximately 217 m resulted in 
the formation of the second and third lagoons and the second 
ice ridge. Numerous blow holes had formed along the second 
ridge by the time the last photo was taken on 17 February. 

The third lagoon was broken up and dispersed overnight and 
a breach had developed in the second ice ridge. During the 
course of the day the second lagoon deteriorated into 
floating brash ice trapped between the fast ice and the 
second ice ridge. There was visible concentration of lag 
deposits along the fast ice located shoreward of the first 
ice ridge, particularly aroirnd the blow hole cones. 

Wave activity during the course of the day broke up the 
second ice ridge and second lagoon. The resulting ice 
blocks and fragments were piled on the lakeward edge of 
the first ice ridge. During breakup a breach in the second 
ridge allowed observation of a wave breaking at a location 
similar to that of the second ridge. 

Little major change occurred in the ice conditions on 20 
February. There was visible fluctuation of a zone of brash 
ice adjacent to the lakeward edge of the fast ice. Blow 
hole cones darkened as the minor snow cover melted from the 
ice surface and sun rot further concentrated lag deposits. 

High air temperatures and easterly winds severely deteri- 
orated the remaining nearshore ice. Sandy lag deposits 
darkened the ice surface and a zone of open water developed 
between the remaining ice and the shoreline. 

Wave activity and a rise in air temperature the previous 
night reduced the remaining ice to grounded blocks approxi- 
mately 1.1 m high and floating brash ice confined between 
the blocks and the shoreline. Storm conditions which 
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23 February 1974 
(CM-14-74-21,22, 
23,24) 



24 February 1974 
(CM-14-74-26,27, 
28,29) 



25 February 1974 
(CM-14-74-30,31, 
33,34) 



26 February 1974 
(CM-14-74-36) 



developed later in the morning forced the remaining ice 
onto the beach and formed lobes of ice along the shore as 
high as 1.7 m. Subsequently the infltix of storm driven ice 
filled these irregularities in the lakeward edge of the 
accreting ice front. The last photo taken on 22 February 
showed a massive, gray white accumulation of shore ice, 
having a relatively smooth lakeward edge in plan that was 
still being augmented by wave tossed ice fragments. 

The wave activity diminished during the night and the first 
slide taken on 23 February showed some minor accretion to 
the fast ice. Due to the incorporation of sediment, the 
ice deposited during the storm had a greyish-brown colora- 
tion in contrast to the lighter colored ice accreted later 
under calmer lake conditions. Photos taken early 23 
February show patches of brash ice adjacent to the lakeward 
edge of the fast ice. Slow accretion throughout the day 
was evident as the white lamination of new ice widened 
along the lakeward edge of the fast ice. 

The first slide taken 24 February showed significant over- 
night ice accretion and a narrow zone of floating slush ice 
just lakeward of the fast ice. Wave tossed ice fragments 
and freezing spray augmented the fast ice, and subsequent 
minor wave activity produced blowholes and an irregular 
shape in plan along the lakeward edge of the fast ice. By 
the time the last photo was exposed 24 February the blow- 
hole cones, located near the edge of the fast ice, were 
nvimerous and well developed. 

The first slide taken 25 February showed approximately 
15 m of overnight ice accretion and slush ice near the 
horizon. The slush ice covered the lake to the horizon by 
early afternoon, and increased in concentration by late 
afternoon. The last slide taken on 25 February showed the 
lake covered with brash ice to the apparent horizon. A 
very dense zone of brash ice was bounded by the fast ice 
and the second breaker zone which was marked by a gentle 
roll on the slushy water surface. 

Due to camera malfunction only two morning slides were 
available for the evaluation of the ice conditions on 26 
February. The lake was covered with ice to the horizon 
and various stages of lake ice development were visible. 
The extinct blow holes were evident lakeward of the greyish 
storm deposited fast ice. Undulations created by freezing 
spray and wave tossed brash ice visible along the ice 
surface marked stationary periods in the lakeward advance 
of the fast ice. A small ridge was forming over the 
position formerly occupied by the second breaker zone. 
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(CM- 15-74) 



Due to camera malfunction the CM-15-74 slide series 
depicting the breakup sequence of the lake ice was lost. 
The first slide taken after camera repair was CM-16-74-1. 
Slides were selected at one week intervals, representing 
the same time of day, until the study area was totally 
ice-free. 



6 March 1974 
(CM-16-74-1) 



The lake was totally devoid of ice to the horizon. The 
remaining shore ice abruptly truncated at the waterline 
forming a nearly vertical lakeward face. The shoreward 
edge of the ice was irregular in plan and covered 
portions of the baseline (seawall) . The entire ice 
body was dark brown in color due to large quantities 
of sediment concentrated on the ice surface by sun rot 
of the ice. A large continuous crack extending parallel 
with the lakeward edge of the deteriorating ice was 
indicative of undercutting by waves and resulted in the 
subsequent collapse of the undercut portion. 



13 March 1974 
(CM-17-74-3) 



Significant narrowing of the shore ice had taken place 
since 6 March 1974. The lakeward edge of the shore ice 
had retreated shoreward exposing the waterline and the 
concentration of sandy lag deposits had further darkened 
the surface of the ice body. 



20 March 1974 
(CM-18-74-5) 



Both edges of the remnant shore ice were very irregular 
in plan. The lakeward edge of the sediment darkened ice 
was located well behind the waterline. 



27 March 1974 
(CM-19-74-5) 



Further visible reduction in the overall mass of the 
deteriorating ice body was apparent in photos taken on 
27 March. The remnant ice now appeared as a series of 
connected ice mounds ranging in height from 0.4 m to 
1.1 m. 



3 April 1974 
(CM-20-74-4) 



By 3 April the shore ice had deteriorated into randomly 
scattered mounds of ice of variable size and shape nearly 
black in color indicating the presence of heavy minerals 
in the lag deposits. 



6 April 1974 
(CM-20-74-16) 



The study area was totally devoid of any ice. 
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APPENDIX B. DESCRIPTION OF SLIDES USED IN ANALYSIS 



Terminology Used in Slide Description 

SLIDE NUMBER: The slide description is presented in order by slide number. 
Each camera monitor slide is coded by roll number, year and exposure 
number. 

TIME: E-AM: early morning 

L-AM: late morning 

NOON : midday 

E-PM: early afternoon 

L-PM: late afternoon 

ACCRETION: A visible lakeward advance of the ice front. 

DETERIORATION: A gradual decay of the nearshore ice complex apparent only after 
viewing several slides in sequence. 

BREAKUP: A rapid decay and dispersal of the nearshore ice usually taking place 
in the period of a few hours. 

ICEFOOT: A small ridge-like ice structure which forms at the waterline. 

ICE LAGOON: A tabular ice formation which shows little vertical relief and is 

confined between the icefoot and first ice ridge or between two ice ridges.. 

ICE RIDGE: A linear ice structure having variable vertical relief and 

usually aligned parallel with the shoreline. A cross-sectional view would 
show a nearly vertical lakeward face and more gently sloping landward face, 

FAST ICE: An ice formation which has greater thickness than an ice lagoon and 
which extends from the icefoot nearly to the location of the first off- 
shore bar. 

SHORE ICE: An ice mass which has been deposited on the beach and is terminated 
near the waterline. 

BRASH ICE: Floating, unconsolidated ice fragments usually formed offshore or 
produced by the break-up of an ice system. 

SLUSH ICE: A floating, unconsolidated ice-water mixture often produced on the 
lake surface during snowfalls. 



1. CM-9-23-74 Extinct blowholes were observed along the first ice ridge, 
(17 January which was discolored due to the concentration of sandy lag 
1974, E-PM) deposits by sun rot of the ice. Grounded, isolated ice 

blocks marked the former position of the second ice ridge. 
The second ice lagoon had been reduced to floating ice 
fragments. 



37 



5. 



CM-10-74-2 
(22 January 
1974, E-AM) 

CM-10-74-3 
(22 January 
1974, L-AM) 

CM-10-74-9 
(23 January 
1974, L-AM) 



CM-10-74-18 
(24 January 
1974, L-PM) 

CM-10-74-20 
(25 January 
1974, E-AM) 



The lakeward edge of the first ice ridge was deteriorated 
sufficiently to expose the first breaker zone. The second 
breaker zone was also defined in the photo. 

Waves breached the first ice ridge showing its mid-ridge 
location relative to the baseline where intersection of 
the ridge with the water surface was evident. 

The first, second and third breaker zones were observable 
lakeward of the ice blocks which were remnants of the 
first ice ridge. Ice fragments derived from the break of 
the first ridge rejuvinated the icefoot. 

The first breaker zone was visible lakeward of the icefoot 
located at the water line. 



Overnight ice formation produced the first ice lagoon , 
composed of floating brash ice, and the first ice ridge. 
A zone of floating slush ice was present lakeward of the 
first ice ridge. 



7. CM-10-74-22 
(25 January 

1974, E-PM) 

8. CM-10-74-23 
(25 Janxiary 

1974, L-PM) 

9. CM-10-74-28 
(26 January 

1974, L-PM) 



The first ice ridge narrowed and its shoreward face steep- 
ened. The first ice lagoon appeared less consolidated than 
in the previous slide. 

The icefoot and first ice lagoon deteriorated while wave 
tossed ice fragments augmented the first ice ridge increas- 
ing its width. 

The icefoot and first ice ridge showed extensive deteriora- 
tion and the first lagoon deteriorated into slush ice. 



10. CM-10-74-29 
(27 January 

1974, E-AM) 

11. CM-11-74-5 
(29 January 

1974, L-PM) 



Overnight wave activity brokeup the first ice ridge and 
rejuvinated the icefoot. The plunge zone, first and second 
breaker zones were observable. 

The first breaker zone and waterline were visible lakeward 
of the deteriorated icefoot. 



12 . CM-11-74-6 
(30 January 

1974, E-AM) 

13. CM-11-74-13 
(31 January 

1974, NOON) 



Due to overnight ice accretion, the lakeward edge of the 
fast ice was located near the first breaker zone. 



Wave activity overnight reduced the fast ice to fragments 
deposited on the beach. The plunge zone and three breaker 
zones were defined. 



14. CM-11-74-30 
(3 February 
1974, L-PM) 



The plunge zone and first and second breaker zones were 
visible lakeward of the deteriorated Icefoot . 
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15. CM-11-74-32 
(4 February 

1974, L-AM) 

16. CM-12-74-12 
(7 February 

1974, L-AM) 



The first breaker zone was defined just lakeward of the 
fast ice which formed overnight. 



New fast ice formed nearly to the position of the first 
breaker zone and blowholes were present along the fast ice 
surface. The first and second breaker zones were visible 
in the slide. 



17 . CM-12-74-18 
(8 February 

1974, NOON) 

18. CM-12-74-20 
(8 February 

1974, L-PM) 



19. CM-12-74-29 
(10 February 
1974, E-PM) 



Slush ice was observed covering the lake surface near the 
apparent horizon. A zone of slush ice had begun to con- 
solidate just lakeward of the static fast ice. 

The accumulation of ice on the lake surface near the 
apparent horizon had consolidated. Brash ice concentrated 
in a band approximately 42 m wide, parallel to the lakeward 
edge of the fast ice. 

Extinct blowholes were noted along the surface of the 
stationary fast ice. The lake was ice covered nearly to 
the apparent horizon and a small ice ridge had formed just 
lakeward of the fast ice. 



20. CM-13-74-13 
(14 February 

1974, NOON) 

21. CM-13-74-19 
(15 February 

1974, E-PM) 



22. CM-13-74-24 
(16 February 

1974, E-PM) 

23. CM-13-74-26 
(17 February 

1974, E-PM) 

24. CM-14-74-2 
(19 February 

1974, NOON) 

25. CM-14-74-3 
(19 February 

1974, E-PM) 



26. CM-14-74-15 
(22 February 
1974, E-AM) 



Breakup produced a large zone of open water bounded on one 
side by the fast ice and on the other by floating ice on 
the lake surface near the apparent horizon. 

The icefoot and fast ice showed further deterioration, and 
sun rot had concentrated sandy lag deposits around the 
extinct blowhole cones. A line of grounded ice blocks were 
observed approximately 216 m from the baseline. 

A small ice ridge was forming just lakeward of the grounded 
fast ice. Numerous randomly scattered ice blocks were 
visible well off shore. 

Overnight ice accretion formed the second ice ridge and 
second ice lagoon. 



The third ice lagoon broke up and the second ice ridge was 
breached by waves. The third breaker zone was visible. 



Ice fragments derived from the breakup of the third ice 
lagoon and breached portions of the second ice ridge were 
forced shoreward by waves. A wave was observed breaking 
through a breach in the second ice ridge. 

Grounded remnant ice blocks marked the position of the 
first ice ridge. 



39 



27. 



28. 



29. 



30. 



31. 



32. 



33. 



34. 



35. 



36. 



37. 



CM-14-74-16 
(22 February 
1974, L-AM) 

CM-14-74-17 
(22 February 
1974, NOON) 



CM-14-74-18 
(22 February 
1974, E-PM) 



CM-14-74-19 
(22 February 
1974, L-PM) 



CM-14-74-21 
(23 February 
1974, L-AM) 



CM-14-74-22 
(23 February 
1974, NOON) 

CM-14-74-23 
(23 February 
1974, E-PM) 

CM-14-74-24 
(23 February 
1974, L-PM) 

CM-14-74-26 
(24 February 
1974, L-AM) 

CM-14-74-27 
(24 February 
1974, NOON) 

CM-14-74-28 
(24 February 
1974, E-PM) 



Storm conditions developed forcing all ice fragments and 
ice blocks onto the beach. The first, second and third 
breaker zones were observed in the slide. 

Large volumes of floating brash ice were forced shoreward 
by NNW winds forming lobes of ice on the shoreline nearly 
2 m high. The entire beach was ice covered and the lake- 
ward edge of the shore ice was very irregular in plan. 

Persisting storm conditions continued to augment the shore 
ice with freezing spray and wave tossed pancake, ball and 
brash ice. The continual influx of ice fragments filled 
the cusps in the. lakeward edge of the shore ice. 

Although storm conditions were beginning to subside, three 
breaker zones were still evident and large quantities of 
floating brash ice continued to augment the fast ice 
complex. 

Three breaker zones were defined as well as slush ice on 
the lake surface near the apparent horizon. The large 
volumes of sediment incorporated in the fast ice deposited 
under storm conditions gave it a grey coloration. The 
lamination of ice which formed along the lakeward edge of 
the fast ice overnight was much whiter making it easily 
distinguishable from the ice deposited during the storm. 

Further accretion of ice was evident along the lakeward 
edge of the fast ice. Two breaker zones were evident. 



The plunge zone and first two breaker zones were observed 
as well as additional ice accretion. 



The plunge zone and first two breaker zones appeared in 
the slide. 



The first breaker zone was evident as well as minor wave 
deterioration of the lakeward edge of the fast ice. 



Further wave deterioration of the ice edge was noted. 
Floating ice fragments were confined between the fast ice 
and the first breaker zone. 

Continued deterioration of the lakeward edge of the fast 
ice resulted in an irregular contour in plan. Further blow- 
hole development indicated the waves had undercut the lake- 
ward boundary of the fast ice. 
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38. CM-14-74-29 
(24 February 
1974, L-PM) 



First and second breaker zones were defined lakeward of 
the fast ice. Brash ice was accumulating between the fast 
ice and the first breaker zone„ Blowhole cones were aligned 
along the outer edge of the fast ice. 



39. CM-14-74-30 
(25 February 
1974, E-AM) 



Overnight ice accretion moved the lakeward edge of the fast 
ice well behind the line of blowhole cones which previously 
marked the furthest lakeward extent of the ice. 



40. CM-14-74-31 
(25 February 
1974, L-AM) 



Little observable change in the variables since the last 
slide. The first and second breaker zones were defined 
and patches of floating ice fragments were seen on the 
lake surface. 



41. CM-14-74-33 
(25 February 
1974, E-PM) 



The lake was covered with unconsolidated slush ice to the 
apparent horizon. The second breaker zone was still 
observable. 



42, CM-14-74-34 
(25 February 
1974, L-PM) 



The lake surface was covered with a dense cover of brash 
ice to the apparent horizon. The second breaker zone was 
defined by a non-breaking roll on the slush ice covered 
lake surface. 



43. CM-14-74-36 
(26 February 
1974, L-AM) 



The lake was ice covered to the apparent horizon. A small 
ice ridge had begun to form over the position of the 
second breaker zone. 



NOTE: The series of slides (CM-15-74) depicting the deterioration and break- 
up of the previous lake ice was lost due to camera malfunction. The 
following slides represented one week intervals beginning with the first 
slide available following repair of the time-lapse photo system. The 
slides were selected to represent as nearly as possible the same time 
of day. 



44. CM-16-74-1 
(6 March 
1974, E-PM) 



The lake was ice-free to the apparent horizon, the only re- 
maining ice was grounded on the beach. The sediment in- 
corporated in the ice at the time of deposition had been 
concentrated on the ice surface by sun rot. The large 
quantity of sediment deposited with the ice under storm 
conditions was evident in the slide. 



45. CM-17-74-3 
(13 March 

1974, E-PM) 

46. CM-18-74-5 
(20 March 

1974, E-PM) 



47. CM-19-74-5 
(27 March 
1974, E-PM) 



Remnant shore ice was confined between the seawall (base- 
line) and the waterline. The edge of the shore ice near 
the waterline had been undercut by waves. 

The remnant shore ice was extremely dark in color due to 
sedimentary lag deposit concentration. Wave activity eroded 
the lakeward edge of the shore ice to a position well behind 
the waterline. 

Further deterioration of the remnant shore ice with increased 
concentration of lag deposits was evident. 
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48- CM-20-74-4 Shore ice had been reduced to isolated patches of ice 

(3 April nearly black in color due to the presence of heavy minerals 

1974, E-PM) in the lag deposits, 

49. CM-20-74-16 No ice was present at the study site. 
(6 April 
1974, E-AM) 
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.APPENDIX C. SOURCES OF POSSIBLE ERRORS IN THE PHOTOGRAPHIC ANALYSIS 



We utilized a high oblique photographic method to monitor the ice con- 
ditions along the shoreline of Lake Michigan. This method, like most, has 
limitations on the accuracy of the numbers that one obtains. This appendix 
describes and where possible quantifies errors that were encountered in the 
analysis. The errors recognized and discussed independently below are: slide 
remounting; marking the location of points on the projected picture; computer 
plotting; approximation of variable positions by a straight line; measurement 
of distances on the plots and position of the apparent horizon in the projected 
picture. 



SLIDE REMOUNTING 

For analysis each slide was removed from the original cardboard mount and 
was remounted on a 2 inch plexiglass square. Alignment of the sprocket holes 
in the color transparency with the guide lines on the plexiglass mount was only 
as accurate as the analyst's eye. Alignment with a given guide line was 
accurate within the width of that guide line (^O.OOl"); hence, this source of 
error was negligible. 

The overall position of the color transparency on the plastic mount was 
also checked. Before the oblique coordinates were assigned to data points on 
the projected image, the position of the photograph's projected principal point 
was verified. If the color transparency had been aligned either too high or 
too low on the plastic mount, the position of the principal point was corrected. 
Typical values for principal point displacement ranged from 0.10" to 0.20" in 
the projected image. 

In order to evaluate the possible error in computed ground coordinates re- 
sulting from a displaced principal point, one slide was analyzed twice. The 
data for the slide were run first without the position of the principal point 
corrected and subsequently with the position correction. The results of the 
comparison indicated that the small errors in principal point location en- 
countered in the study produced differences in the computed ground coordinates 
not exceeding 0.03 m. 



MARKING POINTS ON THE PROJECTED PICTURE 

The computer program utilizes the MSWL as the datum in its distance compu- 
tations. Accurate calculations of ground distances required that the points 
denoting the positions of variables in the oblique image be marked where their 
intersection with the water surface was evident. In this study camera angles 
with respect to the true horizontal were very small and slight variation in the 
camera tilt angle would produce a large change in horizontal distance. The 
small camera depression angle introduced large displacement errors if points 
locating variable positions were marked on the projected photo above the appar- 
ent MSWL. The magnitude of the error in the computed ground distance to the 
variable was a function of the height of the variable above the MSWL and the 
distance of the variable from the camera, both of which determine the tilt and 
depression angles. 
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The relief error (Ax) was calculated for distances from the camera 
ranging from 60 m to 300 m and object heights up to 2 m. These values were 
decided upon as limits since heights calculated in the program never exceeded 
2 m and ice structures were never nearer the camera than 60 m (in the direction 
of the camera' s line-of-sight) . Variables located near 300 m were isolated ice 
blocks and their intersection with the MSWL was easily observable in the pro- 
jected photo thus reducing error in determining their locations. 

The calculated relief correction values (Ax) have been displayed graphically 
allowing interpolation and rapid determination of variable displacement due to 
elevation above the datum. 

We originally intended to calculate the possible displacement errors for 
each month's MSWL separately using the MSWL for the specific month. Figure 14 
was based on the MSWL for February 1974. Calculations based on the MSWL for 
January 1974 showed deviations from the February values were limited to the 
second decimal place for the majority of cases, hence the February values were 
considered representative. The MSWL for the month was calculated by taking the 
average of the MSWL values at both Holland, Michigan and Calumet Harbor, 
Illinois. 

The displacement error due to relief (Ax) is in the lakeward direction 
when ground distance values have been determined from oblique coordinates 
located above the MSWL datum. To avoid large relief errors, ice ridge positions 
were determined at the time of breakup when the intersection of the ice ridge 
with the water surface was evident. 



COMPUTER PLOTTING 

The printer plot subroutines called in the main program constructed two 
cartesian coordinate grids each with a different scale. One grid was scaled 
X = y = 1 inch : 150 ft, and the one used in the analysis. Figure 13, was 
scaled x = 1 inch : 100 ft and y = 1 inch : 150 ft. 

On the plot used here 10 characters were plotted per inch representing 
100 ft along the abscissa resulting in a plotting error of ±5 ft. Along the 
ordinate 24 plotted characters represented 600 ft, hence the maximum plotting 
error in the y direction was ±12.5 ft. Distance measurements made along the 
normal to the baseline had a possible error of ±SQRT (x^ + y^ ) = 
±SQRT 5.02 + 12.52) = ±13.5 ft. ^^'^"^ ^^^^'^ 



APPROXIMATION OF VARIABLE POSITIONS BY STRAIGHT LINES 

The variables were located on the plots in linear arrays that displayed a 
certain degree of scatter about the constructed line representing the average 
position. The points that located the position of the seawall showed the 
least scatter whereas those which located the breaker zones varied the furthest 
from the average line. The points that showed the position of the shore or the 
ice ridges were usually scattered less than those of the breaker zones. The 
scatter of the ground coordinate points from the constructed straight line had 



44 




_J 

CD 
< 

< 
> 



60 120 180 240 300 
DISTANCE FROM CAMERA (METERS) 



FIG. 14. Correction for relief error. 
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Its origin from three sources, the variation incurred in plotting, the natural 
irregularity in the structure itself and the location of the constructed line 
of approximation. 

It was difficult to assign a value to the error incurred by the straight 
line approximation without resorting to statistical methods. Any error that 
did exist was consistent throughout the entire study since all analyses were 
performed by one person. His reasoning for the line placement did not vary 
between cases. We suspect that a variation of 10 to 20 ft was present in some 
cases, but this variation applied to only one or two points along a line whose 
position was represented by six or more points. This was particularly true of 
the breaker zones where points were chosen only where the breaker was evident 
in the picture. As one might suspect, breakers are rarely continuous, and 
breaker zones are marked by a series of breaking waves. 

A regression analysis was later performed on several data sets and the 
correlation coefficients ranged from 0.9 to 1.0 indicating the arrays of ground 
coordinate points (Fig. 13) could be accurately approximated with a straight 
line. 



MEASUREMENT OF DISTANCES ON COMPUTER PLOTS 

Appendix D explains how the distances from the baseline to the lines 
approximating variable locations were measured. The x and y components of the 
baseline normal were determined for each distance measurement . All measure- 
ments of distance on the computer plots were made with a dial caliper accurate 
to 0.001". For the plotting grid used in the analysis the maximum error in 

measurement in the y direction was 0.001 inch * 15 D feet ^^ 0.150 ft. In the 

•^ xnch 

X direction the maximum error was 0.001 inch * 100 fget qj. q.IOO ft. 

inch 

ERROR RESULTING FROM OBSCURED HORIZON 

Since slide CM-10-74-18 had a particularly clear horizon, it was selected 
as the control for evaluation of error in computed ground coordinate values 
resulting from an obscured horizon in the oblique photograph. The distance 
(PTSCRN) measured on the screen between the principal point of the projected 
oblique photograph and the image of the apparent horizon Is used by ICESTUDY 
in the preliminary calculations for the equivalent vertical photograph. Any 
discrepancies in this measurement resulting from an obscured horizon will 
introduce error in the values computed for the principal plane diagram 
(Fig. 2) and subsequent ground coordinate values. 

Increments of ±0.10" were added and subtracted from the control value for 
PTSCRN used in the analysis of slide CM-10-74-18. These plus and minus incre- 
ments represented the possible errors in selection of the location of the 
apparent horizon in a projected oblique photograph with an obscured horizon. 
We are confident that errors in selecting the position of the apparent horizon 
in the oblique representation never exceeded 0.15 Inches. This corresponds to 
a negligible error in the x ground coordinate of less than one foot and less 
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than 8 ft in the y ground coordinate direction. Only 4% to 5% of the slides 
analjrzed could conceivably involve errors of this magnitude. When this 
condition was encountered the selected horizon location was compared with 
an average value for PTSCRN as well as with that for both the preceding and 
following slides to minimize this source of error. 
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APPENDIX D. METHOD USED IN THE DETERMINATION OF MONITORED VARIABLE POSITIONS 

The discussion presented here is the method used in determining the 
location of the monitored variable positions. Figure 15 is used as a basis for 
reference in the discussion. 

1. Construct the baseline (A) through the asterisks which represent the 
location of the seawall. The asterisk located at (57.0, 1247.0) represents the 
position of the north rangepole. 

2. The location of other variable positions on the plot may be approxi- 
mated by lines B, C, D drawn through the plotted characters locating each 
variable. 

(0) = runups (Line B) 

(•) = ice ridge or boundary (Line C) 

(X) = breaker zones (Line D) 

3. Construct the normal (EH) to the baseline (Line A) at point (E) . 

4. Points F, G, H are formed by the intersection of the baseline normal 
and the lines B, C, and D locating the positions of the other variables on the 
plot. From these points of intersection (F, G, H) construct lines HK, GJ, and 
FI parallel to the ordinate of sufficient length to intersect a line drawn 
through point E that is parallel to the abscissa, in this case the horizontal 
printer plot grid line. 

5. Determine the lengths of the X and Y components of the resulting right 
triangle. Ground distances EF, EG, and EH, measured from the baseline, can be 
computed by application of the plotting scale and the Pythagorean Theorem. 
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APPENDIX E. COMPUTER PROGRAM ICESTUDY 



C The ORJECTIVF of this PRORRAM is to OFTFRMINF ground CnORDINAlFS 

C OF ICF RIDGFS. RRFAKFR ZONFS &Mn RR6SH ICF ZONESt THF HEIGHTS OF 

C TCF STRUCTURES AND THE ARFAS OF ICF ACCUMULATIONS FROM HIGH 

C ORLIOUF PHOTOGRAPHS PROJECTED OM A SCREEN THROUGH A SLIDE 

C PROJECTOR. 

C 

C THF VARIABLES USED IN THF PROGRAM: 

C 

C TaIR=TEMPERATURE of THF AIR 

C TH20=TFMPERATURE OF THE SUREACF WATER 

C TAIR, TH20 ARE READ IN DEGREES F AND CONVERTED INTO DEGREES C 

C STICE = ICE CONDITION INDEX - B-UP=BRE AKUP , DETR = DETER lOR ATI ON 

C STAT = STATIC, ACCR = ACCR ET I ON. NONE = NO ICE PRESENT 

C SLTB=PFRPENDICULAR DISTANCE THROUGH THE PRINCIPAL POINT BETWEEN 

C A HORIZONTAL LINE AT THE TOP OF THE PLASTIC MOUNT AND A 

C HORIZONTAL LIME AT THE BOTTOM OF THE PLASTIC MOUNT 

C SCTB=PERPEMDICULAR DISTANCE MEASURED ON SCREEN BETWEEN 

C THF IMAGES OF THE SAME LINES USED IN CALCULATING SLTB 

C SLIDE=ONE HALE THE DISTANCE MEASURED ON THE SLIDE (0.5*SLTB) 

C SCREEN=ONE HALE THE DISTANCE MEASURED ON THE SCREEN (0.5*SCTB) 

C FPROJ=FOCAL LENGTH OF THE SLIDE PROJECTOR 

C FCAM=FOCAL LENGTH OF THE CAMFRA LENS = 50MM OR 1.97 INCHES 

C PTSCRN=DISTANCE PT MEASURED ON THE SCREEN. DISTANCE FROM PRINCIPAL 

C POINT TO APPARENT HORIZON. 

C H=HFIGHT ABOVE THE MEAN STILL WATER LEVEL. 

C SWL=MEAN STILL WATER LEVEL 

C ILFMS=DI STANCE FROM LENS TO SLIDE IN THE PROJECTOR, CALCULATED 

C USING SIMILAR TRIANGLES AND THE GAUSSIAN FORM OF THE THIN LENS 

C FOUATION 

C MnRlZN= HORIZON INDEX, l=EXCeLLENT, 2=GnUl), 3 = EAIK, 4=P(J0R 

C DIP=OIP ANGLE. ANGLE FROM TRUE HORIZON TO APPARENT HORIZON. 

C THFTA=DFPRESSION ANGLE. 

C T=TILT ANGLE OR QO-THFTA 

C PT.PI .TM,IN,TI=DISTANCES ON THE PRINCIPAL PLANE DIAGRAM. 

C RFI = ECAM/H 

C 0=DISTaNCE from PROJECTOR LENS TO SCREEN. 

C GI=V/RFI 

C XBUOY ( ) ,YBUOY( )=COORDINATES OF THE REFERENCE POINTS COMMON TO ALL 

C SLIDES 

C XBZHNEf ) .YRZONE( )=COORDINATES OF THE BREAKER ZONES 

C XRtiNnP( ) ,YRUNUP( )=COnRDINATES OF THE RUNUPS 

C XRI nGF( ) ,YRIDGF( )=COORniNATES OF THE ICE RIDGFS 

C XROYGD( ) ,YBOYGD( )=REAL GROUND COORDINATES OF THE REFERENCE POINTS 

C COMMON TO ALL SLIDES 

C XRZGD( ) ,YBZGD=REAL GROUND COORDINATES OF THE BREAKER ZONES 

C XRUNGD( ) .YRUNGO( )=RFAL GROUND COORDINATES OF THE RUNUPS 

C XRIDGD( ) ,YRIDGD( )=REAL GROUND COORDINATES OF THE ICE RIDGES 

C XICE( ) ,YICE( ) = CO0RDINATES DEFINING THE AREAS OF ICE ACCUMULATION 

C XICtGO( ),YICEGD( )=REAL GROUND COORDINATES OF THE OUTLINE OF THE 

C ICE ACCUMULATIONS 

C ARFA()=ARFA OF THE ICE ACCUMULATIONS AT TIMDAY CALCULATED BY THE 

C COOROINATE METHOD UNITS=SOUARE FEET 

C ARFAM=ARFA IN SOUARF METERS 

C AX,AY=POINT PICKED OFF PHOTOGRAPH AT PEAK OF OBJECT WHOSE HEIGHT 

C IS BEING CALCULATED-USED IN CALCULATING THE EOUIVALENT VERTICAL 

C PHOTO COORDINATES OF POINT AY 

C BX,RY=POINT DIRECTLY BENEATH POINT AX, AY 

C T1=TANGENT OF ANGLE THFT A=TANGENT OF ANGLE TO POINT AX,AV 

C LE=DA=DISTANCE FROM CAMERA TO AX, AY 
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c CB=nisTANce from camera td bx.bx 

C TIMDAY=THe TIME OF DAY IN CATERnRIES E-AM, L-AM, NOON, E-PM, L-PM 

C n IS D PRIMED IN THE VERTICAL HEIGHT FORMULA 

C R IS R PRIMED IN THE VERTICAL HEIGHT FORMULA 

C TS IS THE DISTANCE TS IN THE VERTICAL HEIGHT FORMUA 

C IE ANY OF THE K'S ARF 0, THERE IS NO INPUT FDR THAT PARAMETER 

C K1=NUMBER OF PHOTO RFFFRENCE POINTS (NUMBER OF XBUOY ( ) , YBUOY ( ) 

C PAIRS INPUTFD PER SLIDE) 

C K2=NUMRER OF XBZONE ( ) . YRZONE ( ) PAIRS 

C K3=NUMBER OF XruNUP( ),YRUNUP( ) PAIRS 

C K4=NnMBFR OF XR I DGE ( ) . YR I DGF ( ) PAIRS 

C K7=NUMRFR OF ICE ACCUMULATION AREAS ANALYZED PER SEQUENCE 

C K« = NIIMRFR OF XICF( ) ,YICF( ) PAIRS PFR ICF ACCUMULATION 

C Kq=NUMRER OF 'RY' POINTS PICKED TO USE IN THE ITERATION 

C K10=NUMRFR OF SLIDES ANALYZED IN THIS RUN 

C K15=NUMRER OF POINTS DEFINING SEAWALL LOCATION 

C Klft=K15+l=LAST X.Y BUOY POINT TO BE ENTERED IN THAT DATA SEGMENT 

C THIS LAST POINT LOCATES THE RANGE POLE 

C H1=HFIGHT OF CAMERA AROVE SEAWALL 

C H2=HEIGHT OF CAMERA ABOVE REFERENCE POINT ON RANGE POLE 

C 

C 

C THE PROGRAM THAT FOLLOWS=OBL I QUE PHOTOGRAPHIC INTERPRETATION 

C 

C THE FOLLOWING DEFINES AND 01 MENS lONALI ZES THE VARIABLES 

C 

INTEGER MONTH, DAY, YEAR, PICTO) 

INTEGER HORIZN,TIMDAY 

REAL ILENS,IN 

DI^'ENf, ION XBI)r!Y(50),YRU0Y{ 50) ,XROYGD{ 50) ,YBOYGD,'50) 

DIMENSION XRZ0MF(50) ,YRZ0NE(50) , XBZGD ( 50 ) , YRZGD ( 50 ) , G2 ( 50 ) 

DIMENSION Xr||NUP( 50),YRUNUP( 50 ) , XRUNGD( 50 ) ,YRUNG0(50) 

01 MEN SI ON XRIDGE( 50) . YR I DGF ( 50 ) , XR IDGO ( 50 ) , YR I DGD( 50) 

DIMENSION X I CE( 50) , Y I C E ( 50 ) , XI C EGD{ 50), YICEGD( 50) 

DIMENSION A(50),AREA(50),AREAM(50) 

LOGICAL*! RCOCO /'*', «X', '0', •. «/ 

DIMENSION IMAGE ( 6200 ), X ( 100 ), Y (100 ) 

INTEGER CNTl/0/,CNT2/0/,CNT3/0/,CNT4/0/ 

DIMENSION XI (2000) ,Y1 (2000) ,X2(2000) ,Y2(2000) ,X3(2000) ,Y3( 2000 ), X4 
1 (2000) . Y4(2000) 

DIMENSION Zl(2000) ,W1(2000> ,Z2(2000) ,W2 ( 2000 ), Z3 ( 2000 ) ,W3( 2000 ), Z4 
1 (2000) ,W4(2000) 
C 

C THE ORJECTIVE OF THIS PART OF THE PROGRAM IS TO DETERMINE DISTANCES 
C ON THE PRINCIPAL PLANE DIAGRAM. 

C 
C INPUT THE NUMBER OF SLIDES ANALYZED IN THIS RUN 

RFAD(5,1000) KIO 

IKOUNT=0 
0200 IKOUNT=IKOUNT+1 
C INPUT DATE, SLIDE NUMBER, AIR AND H20 TEMPERATURES AND ICE COND. 

READ ( 5, 901) MONTH, DAY, YEAR, PI CT(1 ) , PICT ( 2 ) , PICT ( 3 ) ,TA IR , TH2n, 
ISTICE 

901 F0RMAT(2I3,I5,4X,3A4,2F10.2,A4) 

C CONVERT TAIR AND TH20 FROM DEGREES F TO DEGREES C 
TAIR=(TAIR-32.0)*0.555A 
TH2 0= ( TH 20-3 ?.0) #0.5556 

WRITE (6,902)MONTH,0AY,YEAR, PICT(1),PICT(2),PICT(3),TAIR,TH20, 
ISTICE 

902 FORMAT ( ' - • , 25X , 31 8, 5X , 3A4, 2F10. 2, 5X, AA ) 
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RFAn (5.900) FCAM 

READ(5,905) SWL.PTSCRN 
C 

C CALCULATE FOCAL LFNGTH OF THE PROJECTOR 

C 

READ (5.906) SLTB.SCTR 

SLIDF = SLTR/2.0 

SCREFM = SCTB/2.0 

n = 9?.n*.SCREEN/(SLIDF + SCREEN) 

ALPHA = 92.0-0 

FPROJ = ALPHA*0/( ALPHA+0) 

RFAn(5,910) HORIZN.TIMDAY 

ILFNS=(n*FPROJ)/(0-FPROJ) 

C=ILEMS/n 
WRITF(ft.800)FCAM,FPROJ.O.SWL.PTSCRN,HORIZN,TIM0AY 

WRITE(ft, 805)1 LFNS.C 
C 

C CALCULATE DISTANCES ON THE PRINCIPAL PLANIF DIAfJRAM 

C 

PT=PTSCRN*C 

\/=PT/FCAM 

H=(62A.546-SWL) 

niP=0.9B*S0RT(H) 

ANGLE=V 

THFTA1=ATAN{ AMGLF )/( 2 .0*3 . 14159/360 .0000 ) 

THFTA=( ( niP/60.0000)+THFTAl )* (2.0*3.14159/360.0000) 

PT=TAN(THFTA )*FCAM 

THETA=THFTA/( 2.0*3.14159/3 60.0000) 

T=(90.0000-THFTA)*(2.0*3.14159/360.0000) 

TT = ri.5*T 

PN=FCAM*TAN(T) 

PI=FCAM*TAN(TT) 

TI=PT+PI 

TM=PN+PT 

1N=TN)-T1 

T=T/( 2.0*3.14159/360.0000) 

WRITE(6,R10)niP.THETAl,THFTA,T,TN,TI,PT,PI,IN 

C 

C THE ORJECTIVF OF THIS PART OF THE PROt^RAM 

C IS To DETERMINE HORIZONTAL OISTANCES 

C 

C 

C IN EACH CASE THE COORDINATES OF THE POINTS ARE FIRST CALCULATED 

C FOR THE EOUIVAI FNT VERTICAL PHOTOGRAPH (COORDINATE AXES ORIGIN 

C AT ISnCFNTER) AND THEN THE REAL GROUND COORDINATES (COORDINATE 

C AXES ORIGIN AT CAMERA) ARE CALCULATED 

C 

C INPUT XBIIOY( ).YBIIOY( ).K1.K15,K16 AND CALCULATE GROUND COORDINATES 

C FOR THE REFERENCE POINTS 



C 



READ ( 5,915)K1,K15,K16 

IF(Kl.FfJ.O) GO TO 100 

H1=3R.63 

H2=15.17 

DO 1 1=1, Kl 

READ(5,920) XBUOY( I I.YBUOYd) 
0001 CONTINUE 

WRITE(6,R20) 
0B20 FnRMAT(' XRUOY(I) YRUOY(I) I Kl ' ) 

DO 2 1=1. Kl 
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WRITE(f>,R25) XBlinY( I ) . YBtlOY { I ) . I ,K 1 
00 0? COMTINIIE 

on ^ 1=1, Ki 

IF ( I.GT.O. ANn.I.I F.K15) H=H1 
IF ( I .GT.KlS.flMn.I .LE.Klft) H=H2 
G1=YBimY( I )*C + PI 

YRiinYi I )=Gi*(Ti/m-Gn ) 

XRtinYl I )=XPI)nY{ I )*c*(Ti/(TI-r,l ) ) 
YROYCnf n=(H*(PI+YRItnY( I ) ) )/FCAM 

xBnYr,n( i )=h*xbiioy( i )/fcam 

0003 CHMTI MtIF 

H= (6?4.546-SWL) 
WRITF(6.R30) 
0R30 FORMAK ' I XBDOYd) YBlJOYd) XBOYGnCl) YBOYGDd)') 

nn 4 1=1, Ki 

WRITE! 6,835)1 ,XBUnY( I),YBIIOY(I ),XRnYGD( I ) ,YRnYr,D( I ) 

0004 CnNTI NIJF 

c 

C IMPDT XBZnNE( ) ,YBZnNF( ) ,K2. AND CALCDLATE THE RROUND COORDINATES 
C FOR THE BREAKER ZONES 
C 
0100 REAn(5.P?5)K2 

IF(K?.EC.O) 0,0 TO 110 

WRITE(ft,R40) 
0840 FORMAT! • I XBZONE(I) YBZONEII) K? •) 

no 11 1=1, K2 

REAni 5.9 30) XRZOME( I ) .YRZONF( I) 

WRITE(A,845)I,XBZ0NE(I),YBZnNF(I),K2 

0011 CONTINDE 

no 12 1 = 1. K2 

G2( I) =YR7.0ME( I )*C + PI 

YBZnNF( n=G2( I )*(TI/(TI-G2( I) ) ) 

XBZnMFC n = XPZONE( I )*C*(TI/(TI-G2 (I ) ) ) 

YR7Gn( I )=(H*( PI+YBZONF( I ) ) )/FCAM 

XRZGDI 1 )=H*XBZONE{ I )/FCAM 

0012 CONTINUE 
WRITE (6, 850) 

0850 F0RMAT(3X,'I XBZONF(I) YBZONEfl) XBZGD( I ) YBZGD(I) G2 ( I ) ' ) 

no 13 1 = 1, K2 
WRITE (6. 8 55) I ,XRZONF( I ) , YBZONE ( I ) . XBZGD ( I ) , VBZGD( I ) ,G2 ( I ) 

0013 CONTINUE 
C 

C INPUT XRtlNUP( ) ,YRIIN1)P( ) ,K3, AND CALCULATE THE GROUND COORDINATES 

C FOR THE RUNUPS 

C 
0110 RFAD(5,935)K3 

IF(K3.F0.0) GO TO 120 

WRITF(6,860) 
ORftO FORMAT( ' I XRUN1IP{I) YRtlMUP(T) K3 ') 

no 20 1=1, K3 

RFAD( 5.940) XRUNUP( I ) ,YRnNI)P( U 

WRITE(6.R65 )I ,XRUN(IP( I ),YRUNUP( I),K3 
0020 CONTINUE 

WRITF(ft,R70) 
0870 FORMATt ' I XRUNUP(I) YRIlNUPd) XRUNGD<I) YRUNRDd ) ') 

no 21 1=1, K3 

G3=YRI)NnP( I )*C + PI 

YRUNUP ( I )=G3*(TI/(TI-G3) ) 

XRilNMP( I )=XRHNIIP( I )*C* ( T I / ( T I-G3 ) » 

YRUNGn( I ) = (H*(PI+YRUNUP( I» ) )/FCAM 
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XRIIMGn( I )=H*XRIINnP ( I )/FCAM 

WR1TE{6.875) I .XRIINIIP( I ),YR()NtlP( 1 ),XRnNr,n( I ) ,YRtlNGD( I ) 

on?i r.nMTiNUF 
c 

C INPUT XRlDr,F( ) ,YRIDr,F( ),K4, ANO C&LCtlLATF THF GRnUND CnnRDINATES 
C FOR THF ICE RinCFS 
C 
0120 RFAn('5,9A5) KA 

IF(KA.FO.O) RO TO 130 
WRITF(f..R90) 
0890 FnRMAT( • I XR I DGF ( I ) . YRIOGF (I ) K4 •) 

nn ^0 i=i,K4 

RFani 5.950)XRinGF( I ) ,YRinGF( I ) 
WRITE( f..885) I ,XRIDGF( n,YRinRF( I).K4 

0030 CnMTIMDF 
WRITE(6.880) 

0880 FnRMAT( ' I XRIOGEd) YRIOGFd) XRIDGOd) YRinGD(I) ') 
nn 31 1=1, K4 
G4 = YRIDGF( I )«C + PI 
YRIOGFI I )=G4*(TI /(TI-G4) ) 
XRinGF( I )=XRinGF( I )*C*(TI/(TI-G4)) 

yRinGn( i ) = (h« (Pi+YRiDr,F( i ) ) )/fcam 

XRinGn( I ) = H*XRinGF( l )/fcam 

WRITE(6.89';)I .XRIDGEf I ),YRIDGF( I ),XRinGn( I ) ,YRIDGD( I ) 

0031 COMTINDF 
0130 COMTINUE 

WRITE (6.4870) 
4H70 FORMAT ( IHl) 

ISnM=l 

KNT1=0 

KNT2=0 

KMT 3=0 

KMT4=0 
C SET UP INFORMATInM RFODIRFD TO CONSTRUCT THF PLOT 

CALL PLOTl (0. IS. 23, 12, 9) 
C PRFPARF GRID AND SET UP INFORMATION RFOUIRED BY PL0T3 TO 
C PLACE POINT IN GRAPH 

CALL PL0T2 (0,300.0,-800.0,3000.0,0.0) 

IF (Kl .EO. 0) GO TO 4876 

no 4871 1=1, Kl 

CNTl = CMTl+1 

KNTl = KNTl+1 

X(I)=XROYGD (I) 

XI (CNTl )=X( I ) 

Zl (KNTl )=X(I ) 

Y{ I )=YROYGD ( I ) 

YKCNTl )=Y( I ) 

WKKNTl )=Y(I ) 
C PLACE PLOTTING CHARACTER IN THF GRAPH FOR EACH (X,Y) 
C COORDINATE POIIMT 

CALL PL0T3(BCn( 1 ) ,X( I ) ,Y( I ) ,1,4) 
4871 CONTINUE 
4876 CONTINUE 

ISUM=I5UM +K1 

IF (K2 .EO. 0) GO TO 4877 

IS=K1+K2 

no 4872 I=ISUM,IS 

CNT2 = CNT2+1 

KNT2=KNT2+1 

X( I )=XRZGD(I-K1 ) 
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4872 
4877 



4873 
4878 



4R74 
487Q 



4891 



X2(CNT?)=X(I ) 

Z2(KNT2)=X(I ) 

Y( I )=YRZGn( I-Kl ) 

Y2(CNT2)=Y( I ) 

W2(KNT2)=Y( I ) 

CALL PLnT3 (BCn(2),X(n.Y(I ).1,4) 

CONTINUF 

CHMTINUE 

ISI)M=1SIIM+K2 

IF ( K3 .PO. n) GO TD 4878 

IS=K1+K2+K3 

nn 4873 I = IStlM,IS 

CNT3 = CNT3+1 

KMT3=KNT3+1 

X( I )=XRIIMGD( I-K1-K2) 

X3(CNT3)=X( I ) 

Z3(KMT3)=X( I ) 

Y{ I )=YRtlNGD< I-K1-K2) 

Y3(CMT3)=Y( I ) 

W3( KNT3)=Y( I ) 

CfiLL PLnT3 (RCn(3) .X(n,Y(I ),1.4) 

COMTIMUF 

CHMTINIIF 

ISUM = ISl)M +K3 

IF (K4 .FO. 0) GO TO 4879 

IS=K1+K2+K3+K4 

nn 4874 I=ISUM.IS 

CMT4 = CNT4+1 

KNT4=KNT4+1 

y ( I >=>'P ! "Gn( I -K 1-K2-K3 ) 

X4(CNT4)=X( I ) 

74(KNT4)=X( I ) 

Y( I ) = Ykin(,f)( I-K1-K2-K3) 

Y4(CNT4)=Y( I ) 

W4(KMT4)=Y( I ) 

CALL PLnT3 (RCn(4),X( I ),Y( n,1.4) 

rOMTI Ml IF 

COMTIMIIF 

PRINT THF rOMPLFTPn GRAPH AMR LARFL 

CALL PLnT4 (20, 'DISTANCE FRHM CAMERA') 

WRITE (6,4891) 

FORMAT ( 1H0,50X, 'DISTANCE LEFT AND RIGHT OF THF CAMERA') 

WRITE (h,5555) 

CALL PLOTl (0,6,23,12,9) 

CALL PLOT 2 ( 0. 600. 0,-lOSO. 0,3000. 0,0,0) 

CALL PLOT 3 (RCD( 1 ),Z1,W1,KNT1,4) 

CALL PLnT3 (RCn( 2) ,Z2,W2,KNT2,4) 

CALL PLnT3 (BCn( 3),Z3,W3,KNT3,4) 

CALL PLnT3 ( RCD(4) ,Z4,W4,KNT4,4) 

CALL PLnT4 (20, 'DISTANCE FROM CAMERA') 

WRITE(6,4891 ) 



THE PURPOSE OF THIS SECTION OF 
THE AREA OF ICE ACCUMULATIONS 



THE PROGRAM IS TO CALCULATE 



INPUT XICE( ),YICF( ),TIMDAY,K7,K8 

RFAn{5,970) K7 
IE(K7.E0.0) GO TO 160 
WRITF(6,730) 
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0730 
1 



OOftl 
740 



0062 



0063 
0161 



FHRMAK 
KR' ) 
DO 60 I 

RFAOCS, 

nn 61 J 

RF/sn( 5, 
WRITF(6 
COMTINl! 
W I T F ( 6 
FORMAK 
• ) 

nn 62 J 

r,6 = YICF 
Y I C F ( J ) 
XICF( J) 
YICFGOI 

xiCEnn{ 

WRITE (6 
CONTI Ml) 
ARFA( 1 ) 
ARFA( 2) 
N3=K8-1 

nn 63 .1 

A( J+l)= 
ARFA( J+ 
COMTIMII 
ARFA( I ) 
8)) 



I 



TIMHAY 



XICF( J) 



YICE( J) 



K7 



= 1,K7 

975)KB,TIMnAY 

= 1,K8 

Q80)XICF( J),YICE( J) 

,735) I .J,TIMnAY,XICF( J),YICF(J),K7,K8 

F 

♦ 740) 

' I J XIC6(J) YICE(J) XICFGD(J) YICEGD<J) 

= ltKR 

( J)«C+PI 

=G6*(TI/(TI-r,6) ) 

= Xir.F( J)*C*(TI/(TI-r,6) ) 

J)=(H*(PI+YICF(J) ) )/FCAM 

J)=H*XICF( J)/FCAM 

,745)1 ,J,XICF(J),YICE( J),XICEr,0( J) ,YICEfiD( J) 

F 

= 0.0 

= 0.0 

= 1.N3 

( (XICFGn( J)+XICEr,D( J + 1 ))*0.50)*( YICFROI J)-YICEGD(J + 1 )) 

1) = ARFA( J) + A( J + 1 ) 

E 

= AREA(K8)-((XICEr,D(l )+XICFGD(KR) ) *0. 5 )« ( YICEGO ( 1 )-YlCEGn(K 



rnNVFRT ARFA( I ) (SOUARE FEET) I MTO SOOARF METERS 

ARFAM(I)=ARFA(n*0.09 29 
WRITF(6,750) 1 ,ARFA( I ) .1 ,AREAM( I ) 
0060 CnMTINME 
160 CnMTINDE 

THIS SECTinN nF THE PRnGRAM CALCULATES HEIGHTS 

INPUT K9 AMn THE PAIRS OF CnnRniMATE PniNTS FOR 
THE HFFGHT CALCIILATinN 



REAn{5,912)K9 
0912 FnRMAT(I4) 

IF(K9.F0.0) GO Tn 170 
nn 4444 II=1,K9 
REAn(5,9R7)AX,AY,RX,BY 
987 FnRMAT (4P10.3) 
G7=PI+C*AY 
YAY=G7*(TI/(TI-G7) ) 
XAX=AX«C*(TI/(TI-G7) ) 
GR=PI+C«RY 
YBY=GR* (TI/(TI-G8 ) ) 
XRX=BX*C«(TI/(TI-G8) ) 
D=G7-Ga 
R=IN+G7 
TS=TI-GR 

HTR=H*( (n«TN)/(R*TS) ) 
WRITE(6,956)G7,G8,D,TS,R.H,TN,HTB 
956 FnRMAT (•-',RF10.4) 

YAYGn=( (H-0.5*HTR)*(PI+YBY) )/FCAM 
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XAXr,n=H*XAX/FCAM 
C CONVERT HEIGHT(HTR) FROM FEFT TO MFTFRS 

HTBM=HTB*0.304R 

WRITE (ft, 5080) 
50H0 FnRMAT(«X'X&Xr,n YAYGO IKOtlNT HTR HTBMM 

WRITE(6,50ai jXAXRD.YAYRn, I KOUNT , HTR, HTBM 
SORl FnRMAT(AX,F10.3.?X,F10.3,5X,I3.SX,F10.3,2X,FlCl.3) 
4444 f, OMTIMtlE 

0170 IF( IKPIINT.LT.KIO) GO TO 200 
0771 COMTINUE 

WRITE (6.5555) 
5555 FORMAT (IHl) 
C 

C THF PLOTS CALLED HERE ARE CUMULATIVE FOR ALL SLIDES ANALYZED 

C IN THIS RUN 

C 

CALL PLOTl (0,6.23.12,9) 

CALL PLnT2 (0, 300. 0,-ROO. 0,3000. 0,0.0) 

CALL PL0T3 (BCn( 1 ) ,X1 ,Y1 ,CNT1,4) 

CALL PL0T3 ( BCD( ? ) , X 2. Y2,CNT2 , A ) 

CALL PL0T3 {RCD( 3) ,X3.Y3,CNT3,A) 

CALL PLOTS (BCn(4),X4.Y4,CNTA,4) 

CALL PLOTA ( 20, 'DISTANCE FROM CAMERA') 

WRITF(6,4891) 

WRITE (6,5555) 

CALL PLOTl (0,6,23,12,9) 

CALL PL0T2( 0,600.0, -1050.0, 3000. 0,0.0) 

CALL PL0T3 (RCn( 1 ),X1,Y1,CNT1,4) 

CALL PL0T3 (BCD(?) ,X2,Y2,CNT2,4) 

CAM PtnT-* (Rrn( ?) ,X3,Y3,CNT3,4) 

CALL PL0T3 {RCD(4) ,X4,Y4,CNT4,4) 

CALL PL0T4 (20, 'DISTANCE FROM CAMERA') 

WRITE(6,4891 ) 
1000 FORMAT( 14) 

900 FORMAT (F10.4) 
0905 FORMAT( 2F10.3) 
906 FORMAT (2F10.3) 
910 FORMAT( I 1 , A4 ) 
091 5 FORMAT ( 314) 
0920 FORMAT(2F10.3) 
0925 FORMAT{ 14) 
0930 FORMAT( 2F10.3) 
09 3 5 F0RMAT(I4) 
0940 FORMAT( 2F10.3) 
0945 FORMAT( 14) 
0950 FORMAT( 2F10.3) 
0955 FORMAT( 14) 
0960 FORMAT( 14) 
0965 FORMAT(2F10.3) 
0970 FORMAT( 14) 

975 FORMAT( 14, A4) 
0980 FORMAT(2F10.3) 
0985 FORMAT( I4,2F10.3) 
0990 F0RMAT(?F10.3) 
0995 FORMAT( F10.3) 
800 FORMAT! •-• ,5F10. 3,15, 4X,A4) 
805 FORMAT(2F10.3) 
0810 FORMAT(9F10.3) 
0825 FOfiMAT(2F10.3,2I4) 

0835 FORMAT( I4,4F10.3) 
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0R46 FHRMAK I4,2F10.3.I4) 

0855 FORMAK IA,5F10.3) 

0R65 FORMAK I4,2F10.3.I4) 

0876 FORMAK I4,4Fl0.3) 

08fi5 FORMAK I4,2F10.3,I4) 

0895 FORMAK I4,4F10.3) 

0705 F0RMAT(2I4,?F10.3,?I4) 

0715 F0RMAT(?I4,4F10.3) 

07?0 FORMAK FIG. 4, 4110) 

0725 F0RMAK4F10.4.2I4) 

735 FORMAK?T5.5X.A4.5X.F10.3,?X,FlO.3.2X, I5,2X,I5) 

745 FORMAT ( 2 I 5. 2X , F 10. 3. 2X , Fl 0. 3, 2 X, Fl 0. 3, ?X. FIO. 3 ) 

7 50 FORMAK 'AREA (' ,13. ')=• ,E20.2,5X, 'ARE AM( ' ,13, •)=' »E20.2) 

0755 F0RMAT(4F10.3) 

0760 FORMAT(2F10.3) 

0770 FORMAK 2F10. 3) 
STOP 
END 
.iFNDFILF 

«sir. 



500 LINES PRINTFD 
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APPENDIX F. SUMMARY TABLE OF CLIMATIC VARIABLES: JANUARY AND FEBRUARY 1974 
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January 1974 
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temperature 
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February 1974 












Air 


temperature 


"C 


Ice 


cond. 


Wati 


er temperature 






-0800 

E-AM 


-1400 

E-PM 


°C 
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1 


-5.0 
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1.7 
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0.6 
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-5.6 


-2.2 


ACCR 


ACCR 


1.1 


* 


1.1 
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-0.6 
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-10.0 


-8.9 
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